Superfluid Heat Conduction and the Cooling of Magnetized Neutron Stars 
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We report on a new mechanism for heat conduction in the neutron star crust. We find that 
collective modes of superfluid neutron matter, called superfluid phonons (sPhs), can influence heat 
conduction in magnetized neutron stars. They can dominate the heat conduction transverse to 
magnetic field when the magnetic field B > 10 G. At density p ~ 10 — 10 1 g/cm the conductivity 
due to sPhs is significantly larger than that due to lattice phonons and is comparable to electron 
conductivity when temperature ~ 10 K. This new mode of heat conduction can limit the surface 
anisotropy in highly magnetized neutron stars. Cooling curves of magnetized neutron stars with 
and without superfluid heat conduction could show observationally discernible differences. 
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Multiwavclcngth observations of thermal emission 
from the neutron star (NS) surface and explosive events 
such as superbursts in accreting NS and giant-flare in 
magnetars provide a real opportunity to probe the NS 
interior (see [lj and [2| for recent reviews). Theoretical 
models of these phenomena clearly underscore the impor- 
tance of heat transport in the neutron star (NS) crust and 
have shown that it directly impacts observations. For ex- 
ample, nearby isolated compact X-ray sources that have 
also been detected in the optical band have a significant 
optical excess relative to the extrapolated X-ray black- 
body emission (a factor 5 to 14) Q. This optical excess 
can arise naturally if heat conduction in the NS crust is 
anisotropic due to the presence of a large magnetic field 
leading to an anisotropic surface temperature distribu- 
tion 0, la, la|- In accreting NSs, the thermal conductivity 
directly affects the observed thermal relaxation time of 
the crust [2j, and the superburst ignition depth and re- 
currence time-scale Q. 

Typically, heat conduction is thought to arise due to 
the flow of relativistic electrons. Electrons are numer- 
ous, degenerate, and provide an efficient mode to trans- 
port heat. Conduction due to lattice phonons has been 
considered, but shown to be unimportant when the crust 
temperature T > 10 6 K [||] and in NSs with low mag- 
netic fields. In this letter we demonstrate that a new 
mechanism for heat transport arising due to the super- 
fluid nature of the inner crust is important at high tem- 
perature and/or large magnetic fields. Here the heat is 
carried by the collective excitations of the neutron super- 
fluid called superfluid phonons (sPhs). 

This mode of conduction is especially relevant in NSs 
with moderate to high magnetic fields (B > 10 13 G). In 
these cases, electron heat transport is very anisotropic 
because electrons can only move freely along magnetic 
field lines. Motion perpendicular to the field is restricted 
because the electron-cyclotron frequency is large com- 
pared to the inverse collision time. In contrast, sPhs, 
being electrically neutral excitations, are not directly af- 



fected by magnetic fields. We find that sPh conduction is 
typically much larger than the electron conduction trans- 
verse to the field and this limits the degree of surface 
temperature anisotropy in magnetized NSs. In NSs with 
low magnetic fields, sPh conduction can still be relevant 
when the temperature is in the range 10 8 — 10 9 K. 

Free neutrons in the inner crust are known to form 
Cooper pairs and become superfluid at temperatures be- 
low the critical temperature T c sw 10 10 K. The superfluid 
ground state spontaneously breaks baryon-number sym- 
metry and gives rise to a new massless Goldstone mode - 
the superfluid phonon (sPh). At long wavelengths, these 
phonons have a linear dispersion relation lj = v s q, where 
v s ~ fcp n /voM, M is the mass of the neutron and k-p-a is 
the neutron Fermi momentum. The thermal conductiv- 
ity of a weakly interacting gas of sPhs can be computed 
from kinetic theory and is given by 



K sPh — - CV V s A s ph 



(1) 



where Cy = 2-k 2 T 3 /(15 v\) is the specific heat of the 
phonon gas and A s ph is the typical mean free path of a 
thermal sPh. We can rewrite Eq.Q] using fiducial values 
of the temperature and the phonon velocity as 
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(2) 
This is to be compared with the thermal conductivity of 
electrons in the inner crust which is approximately 10 18 
erg cm -1 s _1 K _1 when T= 10 8 K [H|. In the presence of 
large magnetic fields electron conduction in the direction 
perpendicular to the field can be several orders of magni- 
tude smaller. Although sPhs are less numerous than the 
electrons they have large mean free paths because they 
interact weakly at long- wavelengths. In what follows we 
will estimate A s ph and its temperature dependence. 

The low energy degrees of freedom in the inner crust 
are lattice phonons (lPhs), electrons, and sPhs. At 



long- wavelengths, the excitations of the ion-lattice with 
mean number density nj, mass number A and charge 
Z, are the longitudinal and transverse lPhs with veloc- 
ities c s = ujp/qTFc and c t ~ 0.7 wp/gez, respectively 
[ll(. Here ujp = \J iire 2 Z 2 ni / (AM) is the ion-plasma 
frequency, qTFo = \/4e 2 /7r ppe is the inverse Thomas- 
Fermi screening length of the electrons, pp c is the elec- 
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tron Fermi momentum, qRz = (4.57r) 
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and 



the inter-ion distance. The electrons are characterized by 
their Fermi momentum pp c = (S^Znj) 1 ' 3 3> m e 3> T 
and form a highly degenerate, relativistic Fermi gas. As 
mentioned earlier, the superfluid phonons (sPhs) are col- 
lective excitations of the neutron superfluid and their 
interactions are weak at low temperature. Their scat- 
tering cross section is parametrically suppressed by the 
factor T 8 //x^ where /x„ = fc Fn /2M is the neutron chem- 
ical potential [12 |. We note that when T ~ T c , sPhs are 
strongly damped as they can easily decay into neutron 
particle-hole excitations. However, since the crust cools 
to T <C T c within several hours of the birth of the NS, 
this processes is exponentially suppressed by the factor 
exp (— 2T C /T) and the mean free path of sPhs is limited 
only by their interaction with lPhs, electrons and impu- 
rities. 

The relevant processes are illustrated in Fig. [TJ (i) 
Rayleigh scattering of phonons due to interactions with 
(compositional) impurities in the solid lattice [131 ]: (ii) 
absorption of sPhs due to their mixing with the longi- 
tudinal lattice phonons which are absorbed efficiently by 
electrons; and (hi) the decay of a sPh into two lattice 
phonons (lPhs). 




FIG. 1: Phonon scattering processes: (A) Rayleigh scattering; 
(B) phonon absorption by electrons; and (C) decay of sPhs 
into lattice phonons. The wavy-line is the sPh, the curly-line 
is the lPh, and the solid line are electrons, respectively. 

At long- wavelength, the cross section for Rayleigh scat- 
tering of phonons by impurities is 
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where q is the phonon momentum and tq is an intrinsic 
strong-interaction length-scale related to the scattering 
length of the neutron-impurity system. This depends on 
the impurity composition and is expected to be about 10 
fm. At low temperature, when q tq -C 1, the mean free 
path of sPhs with thermal energy u>th = 3T is 
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where ni m = 3/(47rd 3 ) is the impurity number density 
and d is the inter-impurity distance. Using fiducial values 
we can rewrite Eq. [4] as 
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where x = d/r a is the diluteness parameter for the im- 
purities and TV is the temperature in 10 7 K. 

We now show that inelastic processes such as sPh ab- 
sorption by electrons and decay to lPhs shown in Fig. [1] 
(B) and (C) are more relevant. To study these processes 
we need a low-energy effective theory which couples the 
sPhs and lPhs. On general grounds, this effective theory 
can only contain derivative terms and at low temperature 
it is sufficient to retain the leading terms. The leading 
order Lagrangian that describes these interactions is of 
the form 
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where </> and £i=i...3 arc the sPh and lPh (canonically 
normalized) fields, respectively. The first term describes 
the mixing between longitudinal lPhs and sPhs, and the 
second and third terms describe the process for the de- 
cay of a sPh into two longitudinal or transverse lPhs, 
respectively. 

The coefficients of this effective theory are determined 
by computing the coupling between neutrons and lattice 
phonons in a manner similar to the calculation of the 
electron-phonon coupling in condensed matter physics 
jl4| . We begin with the fundamental short-range inter- 
action between neutrons and ions given by 



W n l 
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where V(x — x') = 27ra n i S 3 (x — x')/M, is the neutron-ion 
low-energy potential and a n i is the neutron-ion scatter- 
ing length. By expanding the ion and neutron density 
fluctuations in terms of their collective modes we obtain 
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respectively. For typical conditions in the NS crust 
g mix ~ 10~ 3 and A ~ 50 MeV. The specific details of 
this calculation and possible quantum corrections to the 
sPh-lPh interaction will be reported elsewhere 15]. Here 
we employ these results to estimate A s ph- 

Mixing between sPhs and lPhs will damp sPh prop- 
agation because lPhs are strongly damped due to their 
interaction with electrons. To leading order in the small 
mixing parameter g m i x , the mean free path of sPhs with 
energy to is given by 
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where a = c s /v s , and Aiph and riph = Aiph/c s are the 
mean free path and lifetime of the lPh, respectively. 
At low temperature, where the Umklapp processes are 
frozen, we can estimate Aiph by assuming that the lPhs 
decay primarily by producing particle-hole excitations 
(normal processes) in the degenerate electron gas. We 
find that 
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where / cp = Zy/m/AM/ct is related to the electron- 
phonon coupling constant [lq . 

When T> T Um = Z 1/3 e 2 w P /3, the Umklapp pro- 
cesses begin to dominate [T7| and in this regime Aiph-u — 
100 at ~ 2000 fm Q. The corresponding sPh mean free 
path is in the range 10~ 3 — 10~ 6 cm. At intermediate 
temperature we employ a simple interpolation formula 
given by 



-\ph — /« ^1P1 
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where /„ = cxp {-T V JT) [fz| . 

Finally we estimate the mean free path due to the de- 
cay process shown in the Fig. 1 (C). When the sPh veloc- 
ity is greater than the lPh velocity this decay is kincmat- 
ically allowed. The amplitude for the sPh of momentum 
q to split into two longitudinal lPh of momentum q\ and 
(?2 is given by 
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The corresponding mean free path for a thermal phonon 
with energy uj = 3 T can be calculated and we find that 
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where A = A/50MeV, f(a) = 1 - 2a 2 /3 + a 4 /5, and 
a = c s /v s . For the temperatures of relevance we find that 
Aabs < Adccay and A abs < A Ray . Consequently, only the 
sPh absorption process with the mean free paths given 
by Eq. [5] is relevant. 

In Fig. [3 we employ the above result (Eq. |SJ) for the 
sPh mean free path and compare the thermal conductiv- 
ity due to electron, lPhs, and sPhs transport in a typical 
NS crust, at temperatures of 10 8 K and 10 7 K. We show 
the electron contribution parallel (e") and perpendicu- 
lar (e 4 -) to the magnetic field lines for B — 10 13 G and 
B = 10 14 G, according to [18|. Heat conduction due to 
lattice vibrations is primarily carried by the transverse 
modes because of their larger number density and we use 
the results from Q. We can see that, above the neu- 
tron drip point, where a significant number of supcrfluid 
neutrons coexist with the electrons and the lattice, heat 
conduction due to sPh (red dot-dashed curve) is relevant, 
and is always significantly more efficient than conduction 
due to lPhs. The sPh conductivity decreases with depth 
and reaches a minimum when the velocity of the sPh be- 
comes equal the longitudinal speed of sound due to the 
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FIG. 2: Thermal conductivity in a NS crust for T = 
10 s K (left) and T = 10 7 K (right). The dot-dashed (red) 
curve shows the contribution due to sPhs. The thick (thin) 
solid lines show the electron conductivity in the direction 
longitudinal (perpendicular) to the magnetic field lines for 
B = 10 13 G (i) and B = 10 14 G (ii). The dashed lines (yel- 
low) is the lPh contribution. 



resonant mixing. Near neutron drip, where v s <C c s sPhs 
can even dominate over electron conductivity along the 
field. In the direction transverse to magnetic field the 
sPh conduction is shown to be relevant over much of the 
inner crust for fields B > 10 14 G. 



10" 



10 c 



10 



I \ \ 

v N \ 
\ , . rx.A 

; — pole I 

: equator (w/o sPh) 1 

equator (with sPh) 1 






2x1 3 


yr 


J- 


10 8 


. . . . 




10 4 yr 




* i - 
/ .1 . 

y i - 








10 5 yr 




/ r 
-' ^ 

i i 




Ml, I 


&t=: " 


^^7 



1 9 3 4 S 

10 10 10 10 10 



t (yr) 



10 12 10 13 .3 10 14 

p(gcm ) 



FIG 3: Cooling curves of magnetized NS with and without 
considering sPhs. Left panel: Ti, at the neutron drip point vs. 
age. Right panel: T vs. density for fixed times. The polar 
temperature (solid lines) is not affected, but the temperatures 
at the equator with sPhs (dash-dotted lines) and without sPhs 
(dashes) differ significantly. 



To asses how this new mode of heat conduction might 
affect observable aspects of NS thermal evolution we have 
performed magnctar cooling simulations with and with- 
out including sPh conduction. Details regarding the 
simulation and the input physics are described in ear- 
lier work [101 ■ The new ingredients here are the sPh 
contribution to the thermal conductivity and the use of 
an updated So neutron supcrfluid gap obtained from 
Quantum Monte Carlo simulations [2fJ]. The cooling 
curves and t emp erature profiles for a NS model (Model 
A from Ref. []j|) with a poloidal, crustal magnetic field 
of strength B = 5 x 10 13 G at the magnetic pole, are 
shown in Fig. [3l Temperatures at the magnetic equator 
are shown as dashed lines for the model without sPhs and 
as dashed-dotted lines for the model with sPhs. The tem- 
perature at the pole is shown by solid lines and there is 
no significant difference between the two cases. The left 
panel shows Tb, the temperatures at the neutron drip 
point (p«3x 10 11 g cm -3 ), as a function of the NS age, 
t. In the right panel we show the temperature profiles in 
the NS crust as a function of density, for three fixed times 
t = 2 X 10 3 , 10 4 , 10 5 yr. The effect of sPhs is important, 
specially in the inner crust when T sa 10 8 K. In the right 
panel one can see that sPhs partially reduce the other- 
wise strong temperature gradient generated along the in- 
ner crust. As a consequence, the temperature anisotropy 
in the inner crust is limited by sPh conduction. We note, 
however, that an anisotropy generated in the outer layers 
at density below neutron drip cannot be suppressed by 
sPhs. 

Finally, we remark that even the phonon conductivity 
can be anisotropic. Large magnetic fields will induce an 
anisotropy in the electron response because of the gap in 
the electron excitation spectrum due to Landau quanti- 



zation. When this gap is much larger than the temper- 
ature, phonons propagating transverse to the magnetic 
field cannot excite single electron-hole states and this is 
likely to greatly enhance both the lPh and sPh conductiv- 
ity transverse to the field. Rotation will create superfluid 
vortices that will preferentially scatter sPhs propagating 
perpendicular to the rotation axis. Both these effects 
warrant further investigation. 

Our estimate for the sPh conductivity, based on Eq. [TJ 
can be improved by replacing the typical mean free path 
by an appropriate thermal average, and by accounting 
for composition and nuclear structure dependence of the 
neutron-nucleus scattering potential for neutron-rich nu- 
clei in the crust. Nonetheless, our study clearly demon- 
strates that sPhs are important for heat conduction and 
are primarily damped through their mixing with the lon- 
gitudinal lattice phonons. The sPh contribution is signif- 
icantly larger than lPhs and the magnitude of the tem- 
perature anisotropy generated at deep crustal layers is 
likely to be limited by superfluid heat conduction. If fu- 
ture observations allow us to establish a clear correlation 
between surface temperature distribution, magnetic field 
orientation and age, wc may be able to probe the super- 
fluid nature of the NS inner crust. 
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